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The charm quark structure function F% and the longitudinal structure function Ff are directly 
sensitive to the gluon content of proton and therefore are crucial in understanding of proton structure 
function, in particular at low momentum transfer Q 2 and low Bjorken x. In the framework of 
perturbative QCD the charm structure function is calculated in the leading order (LO) and the 
proton structure function is investigated in the next leading order (NLO) at small x region. The 
valence quark distribution is obtained from the relativistic quark-exchange model. The calculated 
F2(x,Q 2 ), F2(x,Q 2 ) and F£(x,Q 2 ), are compared with the present available experimental data. 



I. VALENCE QUARK DISTRIBUTION 

In the evaluation of valence quark distribution, 
we assumed that the nucleon is composed of three 
valence-quarks in the following way [l|: 



peated indices as well as integration over k. <?t (TV" ) 
are the creation operators for quarks (nucleons) and 
■A/?f 1 „ 2 „ 3 are the totally antisymmetric nucleon wave 
functions, i.e. 
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where a designate the nucleon states {P, Ms, Mt} 
and fi stand for the quark states {fc, m s , m t , c}. With 
the convention that there is a summation on the re- 
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The D((Xi, on) depend on the Clebsch-Gordon 

coefficients C^jl and the color factor e ClC2C3 , 
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The 4>{k\, k 2 , /c3, P) are the nucleon wave functions in 
terms of quarks and we write it in a Gaussian form (b 
~ nucleon radius) : 

- - - - /36 4 \' 
d>(ki,k 2 ,k 3 ,P) = 1 — 1 \ 
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We can define the nucleus state based on nucleon 
creation operators, i.e. 

\Ai = 3) = (3!)-^ x QlQ2Q W Ql V Q2 V Q3t |0), (5) 

where ^ Q i Q 2c«3 are ^he complete antisymmetric nuclear 
wave functions (they are taken from Faddeev calcula- 
tions with Reid soft core potential. 

The quark momentum distributions with fixed 
flavour in a three nucleon system are defined as, 



Pp,(k;Ai) 



(Aj = 3]4g4A = 3) 
(At = 3\At = 3) 



(6) 



where the sign bar means no summation on m t and 
integration over k in the fj, indices. By using the above 
definition, we can calculate the quark momentum dis- 
tribution for each flavour. In the above equation we 
use, xi x i Vi cosO), the Fourier transform of the nucleus 
wave function. 

By considering the relativistic correction, the va- 
lence parton distribution at each Q 2 can be related to 
momentum distribution for each flavour according to 
the following equation, 
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Evaluating the angular integration, we find that, 
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with 
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with (yAi) = 1 + jj an d e being the average removal 
energy of the nucleon. A typical ansatz for the parton 
distribution is usually parameterized as, 



where m (M) is the quark (nucleon) mass , fc+ is the 
light-cone momentum of initial quark and eo is the 
quark binding energy. The valence quark distribution 
of a bound nucleon can be derived from the free nu- 
cleon valence quark distribution function by using the 
convolution approximation, 

q v (x,Q 2 ;Ai) = J2 I q v (—,Q 2 ;N)f N/Ai (y Ai )dy Ai , 

(10) 

where f N / A . (y Ai ) are the nucleon momentum distri- 
butions in the nucleus. By taking into account the 
fact that fN/AiiVAi) are large only around , ^ , we 
can write [23] 



'(VAi)' 



,N)^Aq v (x,Q 2 ;A l ) (11) 



x V{x, Q 2 ) = A v rtpx av (1 - x) bv (1 + ^ v x + q v x^) 
where A-p is the normalization factor. 



II. NLO EVOLUTION PROCEDURE 



It is appropriate to use the Mellin and inverse Mellin 
transformation to calculate the NLO parton distribu- 
tions in the (x, Q 2 )-plane [l|. 

In the NLO, F%(x,Q 2 ) is related to the quark, an- 
tiquark and gluon distributions, as follows: 



F£(x;Q 2 ) = x J2 e2{g(z,Q 2 )+g(z,Q 2 )}+^^ 



where <g> means the convolution and it is defined as, where a = 1 + and (rif is the number of active 

flavours) 
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The charm quark contribution to the proton struc- 
ture function F%(x,Q 2 ,m 2 ) has the following form in 

the LO limit, if \ > 1 + (2^) 2 , (note that for small with (A— = 200MeV) 
x, with this condition Q 2 can become less than m 2 , 



4tt/?i lnln(Q 2 /A 2 ) 

/3 ln(Q 2 /A 2 ) ~ PI ln(Q 2 /A 2 ) 
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x g(y, (j,' ) (15) The longitudinal SF,F|>; Q 2 ) is obtained as: 
I 
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Fl(x;Q 2 )=x^4r^ Y, lC q , L ®(q(x,Q 2 ) + q N (x,Q 2 ) + 2C g . L ®G(x,Q 2 )} + Ft(x,Q 2 ,m 2 c ) (18) 
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We use the LO expression for F£(x,Q ,m~), which is 



C q , L = -z, C g , L = 2z(z-l) 



(19) 



IPM-LHP06-19May 



IPM School & Conference on Lepton & Hadron Physics, Tehran-IRAN, May 15 -20, 2006 



3 



the same as equation (3) i.e. 



F c L {x,Q 2 ,m 2 c ) = 2xe 2 c 
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Here we assume the SU(3) flavour-symmetric sea 
quark distributions q = u = d = s=s. In addition 
we consider the sea quark and gluon contributions to 
vanish in the static point [Iq <C Q 2 , ( we set //q = 
0.32GeU 2 ),i.e, 



G(x,n 2 )=0 



q(x,n 2 ) = 0. (21) 



Note that, this is the scale where the above initial 
condition is approximately satisfied i.e. the fig scale 
is determined from the intermediate Q\ scale by evolv- 
ing downward the second moment of the valence quark 
distributions such that the gluon and sea quark dis- 
tributions are approximately zero at the scale. 



III. RESULTS AND DISCUSSIONS 

In figure 1, we present the charm quark contribution 
to the Structure function of proton at ( Q 2 = 6GeV 2 ) 
with different charm mass values i.e. m c = l.lGeV 
(small-dash curve), m c — l.2GeV (dash-small-dash 
curve), m c — 1.3GeV( full curve) and = lAGeV 
(dash curve). The data are those of HI 2] and ZEUS 
collaborations [3j, i.e. the squares (ZEUS, 1997), cir- 
cles (Hl,2000), diamonds (ZEUS,2000) and triangles 
(ZEUS, 2004). Our results are in very good agree- 
ment with the present available data. By reducing 
the charm mass, the charm structure function of pro- 
ton increases but it still passes though the data. It 
also becomes zero for x > 0.1. Obviously the struc- 
ture function becomes larger as we increase Q 2 . The 
calculated F|(x,Q 2 ) by using the gluon distribution 
of GRV (heavy-full curves) with m c = 1.3 shows less 
charm quark contribution to the structure function of 
proton. 

Figures 2 and 3 show the comparison of Q 2 depen- 
dence of charm contribution to the structure function 
of proton for various x values with the available ZEUS 
data (the triangles (1997), circles (2000) and squares 
(2004)). We get a reasonable result with respect to 
the data [2]. 

The Q 2 dependence of longitudinal structure func- 
tion of proton are given in figure 4. The data are from 
HI collaboration experiments: the circles (HI, 2001) 
and triangles (HI, 1996). Again there are a good agree- 
ment between our calculated results and the experi- 
mental predictions data. These show that our gluon 
distribution can reasonably well treat the PGF pro- 
cess. We have observed a linear Q 2 -dependence in the 
unpolarizcd structure function which is in a very good 
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FIG. 1: The charm contribution to the proton structure 
function at (a) Q 2 = 6GeV 2 and (b) Q 2 = WGeV 2 for 
various charm mass. The data are from various ZEUS and 
HI collaborations experiments. 
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FIG. 2: The Q 2 (GeV 2 ) dependence of charm contribution 
to SF of proton for the various x values. The data are 
from ZEUS [2] and HI [3] collaboration experiments. 
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FIG. 3: As figure. 2. 
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FIG. 4: The Q 2 (GeV 2 ) dependence of longitudinal SF of 
Proton. The data are from HI collaboration experiments. 



agreement with the available data. A smooth behav- 
ior was found for the sea-quark and gluon distribu- 
tion functions. The agreement between our uniquely 
gluon and sea quarks distributions and the available 
data are encouraging. In particular we fully repro- 
duced the gluon distribution. We note that the above 
assumption where the gluon and the sea-quark dis- 
tributions can be generated entirely radiatively from 
valence quark may not be valid, especially at small x. 



In summary, by using our recent complete NLO cal- 
culation in the conventional MS factorization scheme, 
we have updated our previous NLO results by includ- 
ing the charm contribution to the SF of proton. We 
have found that the proton structure function has ap- 
proximately the same Q 2 -dependence as the data and 
the whole results are consistent with the available ex- 
periments. Our calculation shows a similar scaling 
violation as the one observed in the experiment for 
the small x. The idea that at low x the scaling vi- 
olation of F^(x,Q 2 ) comes from the gluon density 
alone and does not depend on the quark density which 
was tested in our previous work, is still valid with a 
good accuracy. So, as before, we may conclude that 
the gluon is the dominant source of the parton in the 
small x region. However, it is well known that the 
theoretical interpretation of SF is complicated at low 
Q 2 . Since, in this region the higher twist contribu- 
tions which are proportional to Q~ 2 and Q~ A are not 
included in the DGLAP equations. 
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